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We consider a macroion confined to a cylindrical cell and neutralized by oppositely charged counterions. 
Exact results are obtained for the two-dimensional version of this problem, in which ion-ion and ion-macroion 
interactions are logarithmic. In particular, the threshold for counterion condensation is found to be the 
same as predicted by mean-field theory. With further increase of the macroion charge, a series of single-ion 
condensation transitions takes place. Our analytical results are expected to be exact in the vicinity of these 
transitions and are in very good agreement with recent Monte-Carlo simulation data. 



Properties of charged polymers in solution are intimately 
related to the distribution of small ions around them. A key 
theoretical model for studying this distribution is that of an 
infinite charged cylinder, immersed in a solution containing 
counterions, and confined to a cylindrical cell of finite size. 
When the cell size increases to infinity only some of the coun- 
terions remain bound at a finite distance from the cylinder. 
The remaining ions escape to infinity, leaving behind a dis- 
tribution of ions that compensates only part of the cylinder's 
charge. Furthermore, below a critical linear charge density 
(or, equivalently, above a critical temperature), all the counte- 
rions escape to infinity: The ion density at any finite distance 
from the charged cylinder is zero. 

The existence of a critical temperature, above which all ions 
decondense is predicted within mean-field (MF) theory |lL|2|]. 
We consider here the equivalent two-dimensional (2d) prob- 
lem where ion-ion interactions, as well as ion-macroion inter- 
actions, are logarithmic. Some thermodynamic properties can 
be evaluated exactly in this case, without resorting to the MF 
approximation. In particular, we find that the decondensation 
temperature is the same as that predicted by MF theory - sug- 
gesting that a similar conclusion might hold in the 3d case, 
where ion correlation effects are expected to be weaker than 
in 2d. The same conclusion was pointed out very recently by 
MC simulations in 2d and in 3d @], in which no deviation 
from the MF decondensation temperature was found numeri- 
cally. 

We begin by briefly discussing the MF theory for a charged 
cylinder of radius a in 3d, confined in a cylindrical cell of 
radius R. The MF electrostatic potential depends only on the 
radial coordinate r and obeys the Poisson-Boltzmann equation 
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in which r, the spatial coordinate, was rescaled by the cylinder 
radius: r = r/a, and (j) is the reduced electrostatic potential, 
in units of the thermal energy kT. We assume that counte- 
rions carry a positive charge e and the cylinder is negatively 
charged, with a linear charge density — <?p. This charge den- 



sity enters Eq. Q via t, = /#p, the so-called Manning param- 
eter 0, and Ib = e 2 /kT is the Bjerrum length. The step func- 
tion 9(r) is equal to unity for 1 < f < R/a and to zero else- 
where, and the boundary condition, §'(r = R/a) = 0, enforces 
charge neutrality. Finally, A, is a rescaled fugacity, which does 
not have any physical consequence since changing its value 
merely shifts the MF solution (j) by a constant. The only di- 
mensionless parameters in the problem are thus £, and R/a. 

By defining u — log(f) = log (r/a) and (p = (j) — 2u, Eq. Q 
becomes 
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for < u < L with boundary conditions 
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where L = log (R/a). Equations (|2ji-(|3} can be interpreted 
as describing an ionic solution confined between two paral- 
lel planar surfaces - one at u = 0, another at u = L, having 
surface charges 
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(using units such that /g = 1.) In this equivalent planar prob- 
lem, the surface at u = L is negatively charged and thus always 
attracts the positively charged counterions. On the other hand, 
the surface at u = may be positively or negatively charged, 
depending on For ^ < 1 ions are repelled from the positively 
charged surface, and escape to infinity as L — > °°; For t, > 1 a 
finite fraction of the ions remain bound, so as to neutralize the 
negatively charged surface at u = 0. 

The mapping from cylindrical geometry to a planar one pro- 
vides an instructive way to understand the behavior of the MF 
solution 1 1], but is valid only on the MF level. On the other 
hand, in the 2d case we show that a similar transformation is 
exact, on the Hamiltonian level. 

We begin with the Hamiltonian H n = 2qq'Y!j = \ log(r,/a) — 
g' 2 Liy;l°g| r ; — Tj\, which describes n point-like ions of 
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charge q' interacting with a central disc of charge q and ra- 
dius a in 2d. By analogy with the 3d model, we assume that 
ions are confined to the radial coordinates a < r < R. Charge 
neutrality requires q/q' = n so that, in contrast to the 3d case, 
the number of ions is finite. The partition function is given 
by Z„ = (l/n\)l\ila d^/o 2 * de,-exp(-p7/) where p = 1/kT. 
We define a Manning parameter ^ = fiqq' so that 
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where charge-neutrality is assumed. For convenience, in the 
following we set p = 1 . 

In the grand-canonical ensemble, the partition function 
can be transformed into a field-theory form (as outlined in 
Ref. H) 
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As g' 2 — » the prefactor inside the exponential tends to infin- 
ity. Hence, MF theory [Eq. Q] becomes exact, for any fixed 
value of cj, in the thermodynamic limit n — > °°. In the follow- 
ing, we analyze the canonical partition function Z„ for finite n, 
characterized by the two parameters £, and n (or, alternatively, 
q and q 1 ). To proceed, we note that 
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where m = log(r,/a), C, — exp [(« — l)cjloga], 
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and 



v(m, 9) = —log [2cosht< — 2cos0] 
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The potential v is linear for \u\ ^> 1 being then equal, approx- 
imately, to — |w|. We note that, since < < 2jc is a com- 
pact coordinate, the correction to this linear potential is short- 
ranged. 

A charge q, evenly smeared over the interval, exerts an 
exactly linear potential: 
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It is thus convenient to interpret the linear term in l|8J as com- 
ing from an interaction of the ions with a smeared charge 
qo at u = and a smeared charge q\ at u — L, which re- 
quires —q'(qo — qi) = t, — 2 + Z,/n. Adding the same con- 
stant to go and q\ does not influence the force exerted on the 
ions, and we are free to chose this additive constant such that 
the system is overall charge-neutral, in the following sense: 
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FIG. 1 : The transformation from cylindrical geometry to a problem 
defined on the (t/,9) strip (schematic representation). 



nq 1 = —(qo + qi). With this requirement there is a unique 
choice of qo and q\\ 
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which bears some resemblance to Eq. @. 

So far, we made an exact transformation of the problem 
from cylindrical coordinates into a problem defined on a strip: 
the coordinate u goes from to L, and the coordinate is 
periodic (see Fig. 1). On the (u,Q) strip, ions interact with 
each other through a potential of the form — q' 2 \ui — u%\, aug- 
mented by a short-range contribution. They also interact with 
two smeared charges, qo at u = and q\ at u = L, and the sys- 
tem is overall charge-neutral. Note that any critical property 
of the system should be exactly captured by the long-range 
linear potential term. 

Since we are interested in the behavior when L — > °°, we 
next introduce an approximation, treating the 2d strip as a one- 
dimensional (Id) domain, with a purely linear ion-ion interac- 
tion. This can be thought of as the result of coarse-graining 
on a scale of order 271. Scaling analysis of the partition func- 
tion shows that in the Id model, the value of some observables 
is the same, when L — > °o, as in the 2d problem - for exam- 
ple, the number of ions between u = and u = (XL, for any 
< a < 1 . Therefore we expect the number of bound ions, 
evaluated in the Id approximation, to be the same as in the 2d 
problem |9]. 

In the Id model, the partition function is Z\ c i = 
(l/»0n"=i Jo Ld *, ex P(-^/) where 
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q(x) is the one-dimensional charge density, including the 
boundary charges at and L, and d 2 \|//dx 2 = — 2q(x). Charge 
neutrality ensures that d\\t /dx = outside the interval [0,L] . To 
evaluate Z[ c /, the n particles can be ordered according to their 
position (canceling the 1 jn\ in Z\d). The derivative dx|//ck is 
then equal to — 2qo between and x\ and decreases in a step- 
wise fashion by 2q' at each ion position x,, so that 

Zi d = [ d*i / 6X2- ■■ ( dx n x (13) 

J Jx\ J x n —\ 

exp [-a xi - oci (x2 - x\ ) a n (L-x n )], 

where 0Q = (qo + iq') 2 . Note that this expression could have 
also been obtained directly by writing the partition function 
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0, with a linear electrostatic potential v(x, 9) » — |x| obtained 
from (|9} in the limit L — > °°. 

It follows from Eq. d!3> that Z^/ = /o °/i ° • • ■ of n {L) is the 
convolution of /o, • • • ,/„, evaluated at m = L, where ,/}(«) = 
exp(— OLju) are defined for m > 0. The Laplace transform of 
Zij(L) is thus 
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so that, performing the inverse Laplace transform, 
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Note that neither oca or q depend on L. 

In the limit L — > oo, is dominated by the term k = k* 
having the smallest a^. When the Manning parameter ^ = 0, 
this dominating term is k* = 0; with increase of fc* increases 
in a stepwise fashion, changing by unity at n threshold values 
(where <Xk = oc*-i), 
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Each one of these discontinuities in k* corresponds to a ther- 
modynamic transition. In the following, we analyze the be- 
havior of several quantities at these transitions. More details 
will be presented in a separate publication. 

To evaluate the contact density n(0), it is sufficient to con- 
sider the distribution function of X\ (an ion at x = is neces- 
sarily the closest to the origin) 



» 1 (x 1 ) = (ao-a^)e- (a °- a ** ) - t| 
We thus find that n(0) = do — <Xk* is equal to 
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where 
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Below the first threshold at Z,\ = 1, n(0) vanishes, whereas 
above this threshold it is finite. Therefore the threshold for ion 
condensation is the same as predicted by MF theory. Note that 
the contact density is continuous at % = This is true also at 
each one of the other transitions However the derivative of 
n(0) with respect to \ is discontinuous. 

In the original, cylindrical problem, our result for «(0) 
translates into an ion concentration p(a) = l/(27ta 2 )n(0). In 
Fig. 2 we compare this result with a 2 p(a), as obtained from 
MC simulation of the full 2d problem. Although we used 
the approximate Id model, the agreement between the ana- 
lytical prediction and simulation is very good. In the limit 
n — » oo, the contact density approaches the MF theory predic- 
tion, p(a) -> - l) 2 /(2%a 2 t,). 
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FIG. 2: The contact density, a 2 p(a), as obtained from Eq. i 1 8i 
(lines), compared with MC simulation results |8] (L = 300): n = 1 
(circles), 2 (crosses), 3 (squares), and 5 (diamonds). When Z, < 1 the 
contact density vanishes. 



An exact sum rule, similar to the contact theorem for the 
planar electric double layer |7J], relates the contact density in 
the 2d strip to the number of bound ions: n(0) = — (qo + 
k*q') 2 [in agreement with Eq. 1181 1. This relation is obtained 
by comparing the pressure across the plane u = to the pres- 
sure acting across a plane u = uq, where uq is chosen to be far 
away from both u = and u = L. Since the sum rule is exact 
in both the 2d problem and the Id approximation, equality in 
the number of bound ions implies that p(a) in the cylindrical 
problem, as calculated from Eq. Jl 8i . is exact. 

To evaluate the density at u > 0, the distribution of all ions 
X2,X3, . . . must be evaluated. We find that k* [Eq. dl9i l is equal 
to the number of bound ions, and that for the m-th bound ion, 
the Laplace transform of the distribution function is 
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The Laplace transform of the total particle density is thus 
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This result provides a particularly simple expression for all 
moments of the single-ion distribution in the cylindrical coor- 
dinates, because the (— £)-th moment, 
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where p(r) = n(u) / (Inr 2 ) is the ion density in the cylindri- 
cal coordinates. Figure 3 shows (a/£)(l/r), obtained from 
Eq. ( 1231 (solid lines). Rescaling with /j = a/E, is used to fa- 
cilitate comparison with MC simulation results from Ref. |6] 
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on its inner and outer surfaces. Equation J24> is similar to 
this expression, but there is a correction (second parentheses, 
third term), whose contribution goes to zero only in the limit 
of large n. We expect Eq. J24> to be the exact leading term in 
the electrostatic energy for L — » °°. A comparison with MC 
simulation data (L = 300) [6] is shown in Fig. 4. 

Finally, close to each one of the transition points E, = c^., OCjfc 
approaches ctk-i and, as seen from Eq. dl5> , both and C£_i 
diverge. Concentrating only on their divergent contribution to 
the free energy, we find that 
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FIG. 3: {a/ty{r- 1 ) as calculated from Eq. 03 for n = 1, 2,3,5, and 
100 (alternating solid and dashed lines). Symbols show MC simu- 
lation results 1 6] for the same quantity (L = 300). The dotted line 
shows the prediction of MF theory. 
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FIG. 4: Leading (large L) term in EjinV) [Eq. 1241 . lines], com- 
pared with MC simulation data from Ref. |6] (L = 300, symbols): 
solid line and circles, n = 1; dotted line and crosses, n = 2; dashed 
line and squares, n = 4. The inset shows E/(nL) for n = 15 (solid 
line), together with the MF prediction E/(nL) = J; -1 for i; > 1 (dot- 
ted line). 



on both sides of the transition, where Eq is the leading term 
in L [Eq. (I24H . The leading divergence in the heat capacity 
dE /dT follows as c^/(c^ — <^k) 2 - Scaling arguments, previ- 
ously presented in Ref. 01, are thus in agreement with the 
analytical result. 

In summary, the counterion condensation problem in 2d is 
treated here analytically, taking ion-ion correlations into ac- 
count. A series of single-ion condensation transitions is found 
with increasing in agreement with recent MC simulations 
|0], the first of these occurring at the MF theory transition, 
i; = 1. A possible experimental realization of this problem 
may be obtained with parallel, rod-like polyelectrolytes. Be- 
ing an analogue of the 3d problem with lower dimensionality, 
the 2d model suggests that the Manning transition temperature 
in 3d is exact even in the presence of ion-ion correlations. 

We acknowledge discussions with A. Naji and R. R Netz, 
and thank them for sharing with us their simulation data prior 
to publication. This research was supported in part by the 
National Science Foundation under Grant No. PHY99-07949. 



(symbols). The agreement is good but not perfect - deviations 
reflect the effect of using the Id model as an approximation to 
the 2d model on the strip. As should be expected, agreement 
is perfect in the MF limit (n — > °°) as well as in the opposite 
limit, n = 1 . 

We next evaluate the electrostatic energy. To leading or- 
der in L, the free energy in the Id model is — logZ^/ = (Xk*L, 
which corresponds to a free energy (<x* k — q\)L in the original 
cylindrical problem. The mean electrostatic energy E is found 
by taking a derivative c^(d/dc^), which is equivalent to p(d/d(3) 
by virtue of Eq. (|5}, and yields 

? 

E = ~(n-k*)(n-k*+l)L. (24) 
n 

One may expect to find E = q' 2 (n —k*) 2 L, the electrostatic 
energy of a cylindrical capacitor having charges ±q'(n — k*) 
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